Tracking fast neutrons 
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Based on elastic collisions, the linear momentum of a fast neutron can be measured from as few as 
two consecutive recoil ion tracks plus the vertex position of the third collision, or 'two and half ion 
tracks. If the time delay between the first two consecutive ion tracks is also measured, the number 
of ion tracks can be reduced to one and a half. The angular and magnitude resolutions are limited 
by ion range straggling to about ten percent. Multi-wire proportional chambers and light-field 
imaging are discussed for fast neutron tracking. Single-charge or single-photon detection sensitivity 
is required in either approach. Light-field imaging is free of charge-diffusion-induced image blur, 
but the limited number of photons available can be a challenge. 1 H, 2 H and 3 He could be used for 
the initial development of fast neutron trackers based on light-field imaging. 



I. INTRODUCTION 

Recently, several experiments were reported that si- 
multaneously measured the directions and energies of fast 
neutrons from recoil ion tracks [1—4]. In the experiments, 
*H or He nuclei were used as the ion targets, the ion- 
ization tracks of MeV recoil protons or a-particles were 
recorded using either electronic or optical readouts. In 
addition to fast neutron counting, the directions and en- 
ergies of fast neutrons were inferred from the ion tracks. 
The collisions were predominantly elastic. A neutron- 
tracking scheme based on 3 He absorption was also de- 
scribed [1]. 

It is well known that a single elastic encounter between 
a neutron and a nucleus is insufficient to transfer all of 
the linear or translational momentum of the neutron to 
the nucleus, except for the rare occasion of a head-on col- 
lision between a neutron and a proton [5-10]. Here we 
show analytically that, based on elastic collisions, the lin- 
ear momentum of a fast neutron can be measured from as 
few as three consecutive elastic encounters between the 
neutron and a low-Z target consisting of 1 H, 2 H or 4 He 
atoms. The effects of an electromagnetic field or grav- 
ity on fast neutron motion can be ignored here. A fast 
neutron trajectory is a line segment in between two con- 
secutive collisions, which can be determined by localizing 
the vertices of the collisions. The use of a higher- Z tar- 
get, such as 12 C, is hindered by low detector sensitivity 
to ions with energies below 100 keV. The energy transfer 
from a neutron decreases when the atomic number A of 
the target increases [6-9]. 

Besides a head-on elastic collision with a proton, a fast 
neutron can also transfer its full linear momentum to a 
target nucleus in a single collision when it is absorbed. 
3 He, 6 Li and 10 B are of interest because of their rela- 
tively small masses [10]. In each case, the corresponding 
compound nuclei 4 He*, 7 Li*, and n B* are not stable and 
disintegrate into two ions that can be tracked, and the 
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fast neutron momentum can be deduced. In the case of 
n B*, 94% of the decay branch also emits a Doppler-shift 
7-ray, which may be used to obtain neutron momentum. 
The neutron absorption cross sections are less than 1 
barn for MeV neutrons, high detection efficiency would 
require a large amount of target mass. Indeed, neutron 
absorption by 3 He, 6 Li or 10 B is widely used for ther- 
mal neutron detection since the absorption cross sections 
scale inversely with the speed of a neutron [10-12]. Al- 
though the present discussion of fast neutron tracking 
primarily focuses on elastic collisions, the measurement 
techniques below should also be applicable to absorption- 
based neutron tracking. 

Below, we first describe fast neutron tracking principle 
in Sec. II. Next, we show that the measurement errors 
are limited by ion range straggling. In Sec. IV, we discuss 
neutron tracking using multi-wire proportional chambers 
and light-field imaging. In the following section, compar- 
isons between elastic collision based methods and absorp- 
tion based methods are made. A brief summary is given 
towards the end. 



II. NEUTRON-TRACKING PRINCIPLE 

A fast neutron (MeV in kinetic energy) moves in a 
straight line and makes a sudden turn occasionally when 
it gets in a close proximity to ~ 10~ 15 m of another 
nucleus. For a low-Z nucleus like 1 H, 2 H, 3 He, 6 Li, or 10 B, 
elastic collision and absorption are most likely during the 
close encounter. The details of nuclear interaction are 
not needed for an elastic collision and the trajectories of 
the scattered neutron and the recoil nucleus (it is used 
with 'recoil ion' interchangeably) can be determined from 
the conservation of momentum and energy within the 
framework of classical mechanics [6-9] . The momentum 
of the scattered neutron, as well as the momentum of 
the recoil nucleus, is random. The angular distributions 
of the scattered momenta do depend on the details of 
nuclear interaction, such as the type of nucleus, nuclear 
spin and neutron energy. 

In the first scenario, which solely relies on the recorded 
positions of recoil ion tracks, three ion tracks are suffi- 
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cient to measure the linear momentum of the fast neu- 
tron. We shall assume that the track lengths can be 
used to derive recoil ion energies. The assumption is 
valid to within the straggling ranges of ion stopping. The 
straggling ranges are within a few percent of the mean 
stopping ranges of MeV ions. In Fig. 1, the measured 
quantities are recoil ion momenta Qi, Q2, and Q3. The 
initial neutron momentum is Pq. The subsequent neu- 
tron momenta are labelled as Pi and P2. Momentum 
conservation gives 

P = Qi + Pi, (1) 

with 

P 1 = Q 2 + P 2 . (2) 




It is clear that the full recording of the third ion track Q3 
is not necessary. The collision vertex C, together with the 
location of B, is suffice to obtain e2. In short, the vectors 
Qi, Q2, ei and e2, or '2 and half ion tracks, are sufficient 
to determine Pq. If one approximates the ion tracks by 
straight lines, the set {Qi, Q2, ei, e 2 } is derivable from 
a five-point measurement, the three vertices A, B and C, 
together with the end points of the first and second recoil 
ion track. 

In the second scenario, the five-point measurement can 
be reduced to a three-point measurement, if the neu- 
tron time-of-flight (TOF) can be measured in between 
ion tracks. Specifically, if the time delay between the 
first (point A in Fig. 1) and the second (point B) col- 
lision, symbolized by n, is measured, then t\ uniquely 
determines the magnitude p\ and the third collision is 
not necessary. Only the vertices of the first and second 
collisions are needed to deduce ei. The set {Qi, ei, n} 
determines Po through Eq. (1). 

One may ask whether consecutive Compton scatter- 
ings of gamma (7) rays can be used for tracking 7 rays. 
The fact that 7 rays always move at the speed of light 
necessitates the five-point measurement. The stopping 
powers for Compton electrons are smaller than those of 
ions in the MeV energy range. Electron tracks can not be 
approximated by straight lines. These differences make 
7-ray tracking unique in its own way and 7-ray tracking 
will not be discussed further. 



FIG. 1. Fast neutron linear momentum measurement based 
on three recoil ion tracks. A, B and C are the three vertices 
of three consecutive elastic collisions. The ion momenta are 
Qi, Q2, and Q3. The initial neutron momentum is Po. The 
neutron momentum after the ith collision is labeled as Pi. 
Qi, Q2, and Q3 are in general not on the same plane. But 
the triplets {Po, Pi, Qi} and {Pi, P2, Q2} are each coplanar 
because of momentum conservation. 

The directions of Pi and P2, symbolized here by the 
unit vectors ei and e2, can be obtained from the locations 
of the vertices A, B and C. The momenta satisfy Pi = 
pi ei and P2 = P2 &i respectively with p\ and P2 being 
their magnitudes. The dot product of ei with Eq. (2) is 

Pi = Q2 • ei +P2ei ■ e 2 . (3) 

The dot product of e 2 with Eq. (2) is 

Pi^i ■ e 2 = Q2 • e 2 +P2- (4) 

Solving Eqs. (3) and (4) for pi and P2, 

Q2 ei (Q 2 ■ e 2 )(ei ■ e 2 ) 



and 



Pi 



P2 



l-(e!-e 2 ) 2 

Q2 e 2 - (Q 2 • ei)(ei ■ eg) 
1 - (ei ■ e 2 ) 2 



(5) 



(6) 



III. MEASUREMENT ERRORS 

The error in a vector measurement can be broken down 
to the error in magnitude and the error in direction, 5K.Q 
= <5k e + k <5e . Here K = fcoe is the vector. <5k e 
describes the error in magnitude and ko^eo the error in 
direction, since <5eo-eo = 0. For example, for a length 
vector based on the measurement of two end points, ro 
= T2 - T±, the relative error in length (e, ) is given by 



ro 



V251 
ro 



(7) 



where we assume that the errors in the three coordinates 
are the same ~ 81 for both ends ri(xi, yi, z\) and Y2(x2, 
2/2, Z2). That is, 8x\ = 8y\ = ■ ■ ■ = Sz2 = 81. The error 
in direction (e^) is given by 



\r 5e \ 



V 6e ■ 6e . 



It is straightforward to show that 

261 

£0 = • 

ro 



(8) 



(9) 



More generally, if the errors in different coordinates 
vary, one may define 



6l 2 



8x1 + 8x1 8yl + 5yl 8z\ + 8z\ 
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then Eqs. (7) and (9) are valid as the upper bounds on 
the errors using the generalized definition of 51. 

The experimental implications are that a.) the ratio 
of an ion track length to the spatial resolution of a track 
recording method determines the errors in magnitude 
and direction of an ion momentum Q; and b.) 51 <C 
is desirable for small errors. 

In the fast neutron tracking problem, fortunately, there 
is a natural length scaling, A 3> L ^> 6L s i ~ 5L st , as 
shown in Fig. 2. Here A is the mean- free-path of neu- 
tron elastic collision. L is the recoil ion range. 5L s i is 
the ion longitudinal straggling range and 5L st the lat- 
eral or 'transverse' ion straggling range. This length 
scaling implies that errors in neutron directions, Pi and 
P2 in Fig. 1, which are proportional to 51 /X (much less 
than 1%), are much smaller than the errors in ion direc- 
tions, Qi and Q2, which are proportional to 51/ L (around 
10%). Here 51 stands for the detector spatial resolution. 
In the analysis below, we will neglect the direction er- 
rors in neutron momenta after the first collision, Fig. 3. 
The error in the neutron momentum Po comes from the 
directional and magnitude uncertainties in Qi and the 
magnitude uncertainty in Pi. 



to its energy through 
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FIG. 2. A natural length hierarchy exists in the fast neu- 
tron tracking problem, A > L > 5L s i ~ 5L s t, with A being 
the mean-free-path of neutron elastic collision, L the recoil 
ion range, SL„i being the longitudinal and 5L„t the lateral 
('transverse') ion straggling ranges respectively. The calcula- 
tions are for 1 H, 2 H, 4 He with a molecular density of 2.69x 
10 20 per cm 3 (or 10 bar pressure at STP). The ratios of the 
lengths are approximately invariant to density changes. The 
neutron mean free paths are based on cross sections given 
by [13]. The ion ranges and range straggling are outputs 
from SRIM calculations [14]. 



5E X 

PlC 



(12) 



with Ei — (71 — l)mic 2 . Here we keep the relativistic for- 
mulation for the ion. /?i is the ratio of the ion speed to the 
speed of light c, and so forth. Since Ex = E\{L{), the en- 
ergy is a function of the ion range, 5Ei = 5L\{dE\/dL\). 
Therefore, 



5qi = 



5 ^ dEx 
fixe dLx 



(13) 




FIG. 3. Because fast neutron mean free path is much greater 
than recoil ion range, a good approximation is that the er- 
rors in the directions of neutron momenta after the first col- 
lision can be neglected. Therefore, the error in the incoming 
neutron momentum Po comes from both the directional and 
magnitude uncertainties in Qi and only the magnitude error 
in Pi. 

Two factors contribute to the ion range uncertainty 
5L X - One is the instrumental error due to the finite de- 
tector spatial resolution, 51. The other is the longitudinal 
straggling. Since the two uncertainties are not correlated, 
one has 



5Lx < ^25l 2 + SL^Qxy 



(14) 



where we use 5L s (Qx) = max [L s /(Qi), L st (Qx)} for the 
upper bound on the range straggling for the first re- 
coil ion momentum Qi. There is no factor of 2 in the 
5L s i(Qi) 2 term because only the end point of the ion 
track is affected by straggling. 

From above, the magnitude error in Qi is 



We begin by examining the error in the first recoil ion 
momentum, <5Qi, which can be determined from two po- 
sitions, the vertex A and the end point of the recoil ion 
track. 



r5Qi = SqJx + qxSfx, 



(11) 



here we use Qi = gifi and fi is the unit vector for the ion 
momentum Qi. The momentum of a recoil ion is related 



,~ \ Sqx 5Lx dEx 

e r (Qx) = = , 

qx 2Ei dLx 



(15) 



which is accurate to the order 0(f3f) ~ 1%. 

The direction error in Qi, which is the same as the 
error in fi , 



£0(Qi) = £0(fi) 



y /46P + 26L a {Q 1 ) 2 



(16) 
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From Eqs. (14) and (16), there is no significant ben- 
efit to achieve an instrumental spatial resolution much 
below the range straggling, 81 <C 5L s (Qi). Matching the 
spatial resolution and ion range straggling 61 ~ <5L s (Qi) 
is sufficient. Fig. 4 shows that, by adjusting the stop- 
ping gas pressure, a detector resolution of 200 /im can be 
matched with the ion straggling range. 




4 6 

He4 ion energy (MeV) 

FIG. 4. Optimal detector pressure is to match the straggling 
range of an ion with the detector spatial resolution, or \f28l < 
<5L s (Qi) as in Eq. (39). 

Next we examine the error in Pi as a function of Q2, 
ei and e2. It is convenient to replace the pair {ei, 62} 
by a transformed pair {h, £ = sin O2} with the new unit 
vector h defined as 



and 



gi - (ei ■ e 2 )e 2 
C 



C 2 = 1 - (ei • e 2 ) 2 



(18) 



02, as defined in Fig. 1, is the angle between the vectors 
Pi (or ei) and P2 (or e2). In the new notation, Eq. (5) 
for pi becomes 



Pi 



Q 2 h 

C ' 



(19) 



We also have h-ei = £ and he 2 =0. In general, the 
variations of <5Q2, 8e% and Se2 are not independent of 
each other since they share a common point B. Due to the 
above length scaling, we use the following approximations 



e^(Pi) = 8ei = 8e 2 = 0. 



(20) 



Therefore, 



< 



8Q 2 ■ h (5Q 2 ■ h 

c c 
c 2 

8q% + q 2 8f 2 ■ 8f 2 



(21) 



As a result, 



£r(Pl) 



g2 
Pit 



£ 2 (Q 2 ) + d(Q 2 ) 



(22) 



here er.(Q 2 ) and e^(Q2) are similar to Eqs. (15) and (16), 
the error in Qi. 

At last, we can calculate the error in Po as a function 
of relative amplitudes and directions of Pi and Qi, that 
is, <5Pq = 5Pi + <5Qi- In terms of the components, 



Sp e +p 8e = Spie-i + <5Qi, 



(23) 



here again, we have neglected the error <5ei. This equa- 
tion can be broken down to the magnitude error 



8pl = Spjie^eo) 2 + (SQ 1 ■ e ) 2 
< ^cos6 2 + 5Qi • 6Q1 
= 8p\ cos Oj + 8q\ + qfSf! ■ , 

and the angular or direction error 

Se ■ 5e a _ 1 ei ■ 5e a ^ SQi ■ 5e 
\Se \ p 1 |£e | \Se \ 

The direction error can be simplihed to 

Se ■ Se < \(Spl sine 2 + 5q\ + ql5h ■ 8h) 
Pa 

From Eq. (24), one has 



(24) 



(25) 



(26) 



£r(Po) < -^/p 2 e r (Pi) 2 cos 2 9i + g 2 e r (Qi) 2 + g 2 e (Qi) ; 

(27) 



(17) From Eq. (26), 



£0(Po) < -VP 2 er(Pi) 2 sin^ei+( ? 2 e r (Qi) 2 + (? 2 e (Qi) 2 
Po v 

The discussion thus far is on the five-point method. 
For the three-point method when the neutron TOF (ti) 
between A and B in Fig. 1 is available, pi, instead of 
Eq.(5) or Eq.(19), is now given by 



Pi 



Lab 



m c 



TlC 



1 



• Lab \2 

' TIC ' 



(29) 



The second fraction is the relativistic correction, which is 
on the order of f3f with 0i being the ratio of the neutron 
velocity to the speed of light. For a 10 MeV neutron, the 
relativistic correction is about 1%. Therefore, with 1% 
precision, it is sufficient to use the approximation 



L 



Pi 



AB 



TlC 



1 ( L AB 

2 V nc 



(30) 



Therefore, 



e r ( 
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As in the five-point method, we neglect the errors on the 
order of SL AB /L AB ~ 51/X < 0(1%), then 



Tl 



(32) 



The component errors in Po are still expressed by 
Eqs. (27) and (28), with the understanding that e r (Pi) 
is now given by Eq. (32). One may also try to use the 
TOF method to measure the ion momenta and obtain 
an error in the form given by Eq. (31). One now notices 
that the error due to first term SL/L in Eq. (31) could 
be comparable or larger than the error due to the time 
measurement. Therefore, due to ion range straggling, 
the TOF measurement does not significantly reduce the 
errors in measurement of Q;. 



A. A neutron track example 

To obtain a quantitative understanding of the errors 
described by Eqs. (27) and (28), we need the details of 
fast neutron scattering. The results are well known, in 
particular, when the angular anisotropy of the scattering 
in the Center-of-Mass (CM) frame, due to nuclear spin 
and other internal degrees of freedom, are neglected [5- 
8, 15]. The ratio of the neutron energies before and after 
the nth scattering is given by 



E„ 



A 2 



1 + 2A cos 9„ 



1 



o 



K-i {A + l) 2 

and the scattering angle is given by 

A cos 6„ + 1 

cos©„ = 



y/A 2 + 2Acos9 n + 1 



(l-cos0 n ), (33) 



(34) 



Here following [5, 9], the a parameter is defined as 
a = 4 Aj (1 + A) 2 . This definition of a is somewhat arbi- 
trary [8]. O n is the counterpart of ©„ in the CM frame. 
While it is tempting to average Eqs. (33) and (34) and use 
them to calculate the errors described by Eqs. (27) and 
(28). We do not use this approach here. A more rigorous 
treatment of angularly averaged errors, like Eq. (27), will 
require integrations over scattering angles Oi and 6 2 , or 
equivalently 9i's in the CM frame [6, 15]. The angular 
average of a term like (p 2 e,-(Pi) 2 cos 2 ©i) in Eq. (27) is 
not the same as the product of the averages of its factors, 
(p 2 )(e r (Pi) 2 )(cos 2 ©i). Individual neutron tracks, how- 
ever, do not depend on ensemble-averaged properties. 
The neutron track example here satisfies 



1 



a 
2' 



(35) 



which has nevertheless the same energy ratio as the av- 
erage derived from Eq. (33). This result is equivalent to 
pick 6 n — 7r/2. Correspondingly, 

cos©„ = , (36) 

v^TT 

which is different from 2 / (3A) when averaging Eq. (34) . 



B. Tracking errors 

We now estimate the measurement errors of the neu- 
tron track that satisfies Eqs. (35) and (36). Eq. (15) 
gives 



-(Qi) 



< 



\/2&L a (Qi) dE 1 

2E~i ~dL~x 
V2 SL s (Qi) L x dEi 

~~2 L~ x E\~dL~x 
V2^MQi) 

2 L x 



(37) 



To obtain Eq. (37), we use the approximation 
(Lx/Ex)(dEx/dLx) < 1, which is valid for Ex above the 
Bragg peak of energy loss dE/dL. Based on SRIM calcu- 
lations for 1 H in Methane gas [14], the Bragg peak is at 
Ex = 0.06 MeV. For 2 H in 2 H gas, the Bragg peak is at 
Ei= 0.11 MeV. For 4 He, the Bragg peak is at Ex = 0.65 
MeV. Eq. (16) gives 



e*(Qi) = 



2SL s (Qx) 
Lx 



Here we have chosen 



V261 < SL s (Qx)- 



(38) 



(39) 



The '<' sign is because two ions are involved. Therefore, 
\/25l ~ mini(f5-L s (Qj)). The results are shown in Fig. 5 
for a variety of conditions. The ratios are close to 4% 
for longitudinal straggling, and less for lateral straggling. 
The results lead to e r (Qi) of about 3% and e«,>(Qi) of 
about 8%. The results also apply to the second recoil 
ion, that is, e r {Q,2) ~ 3% and e^(Q2) ~ 8%. The error 
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FIG. 5. Ratio of the range straggling to the total range at dif- 
ferent ion energies. The solid lines are for longitudinal ranges. 
The dotted lines for lateral ranges. 1 H is in methane gas. 2 H 
is in deuterium gas. The results are not sensitive to pressure. 



in Pi, as given by Eq. (22), is 

6r(Pl) 



yz(A 2 + 1) 



A + l 



^(Q 2 ) + d(Q 2 ) 



(40) 



6 



Therefore, for 1 H recoils, e r (Pi) = 8.5%. For 2 H recoils, 
e r (Px) = 9%. For 4 H recoils, e r (Pi) = 10%. The errors 
in P , as given by Eqs. (27) and (28), are 



£r(P ) < 



V2A 2 + 1 
A+l 



e 2 (Q 2 )+d(Q 2 ), (41) 



and 



e*(Po) < 



V3A 
~A+1 



e?(Q 2 ) + d(Qa). 



(42) 



For X H recoils, e r (P ) = e (P o ) = 7.4%. For 2 H recoils, 
e r (P ) = 8.5%, e (P o ) = 10%. For 4 He recoils, e r (P ) 
= 10% and e (P o ) = 12%. 

If TOF is used, e r (Pi) is free of errors from ion strag- 
gling. If e r (Pi) is small compared with ion straggling, 
then 



e r (P ) =e (P o ) < 



V2A 

a + i 



e 2(Q 2 )+e 2 ,(Q2). (43) 



For X H recoils, e r (P ) = e^Pn) = 6%. For 2 H recoils, 
e r (P ) = e (P o ) = 8%. For 4 He recoils, e r (P ) = e (P o ) 
= 10%. These results are not significantly different from 
above. Therefore, a TOF accuracy 



5n 



< 



Tl 



e?(Q 2 



•4(Q 2 ) 



(44) 



is sufficient. 



C. Tracker dynamic range and size 

The spatial dynamic range (2)) of a fast-neutron 
tracker may be defined as 



(45) 



the ratio of the neutron mean free path to the spatial res- 
olution of the track recording instrument. From above, 
SI = 5L s /^/2, the spatial resolution is chosen to be com- 
parable to the ion range straggling. The spatial dynamic 
range is plotted as a function of ion energy in Fig. 6 for 
three different ions. £> is in the range of 10 4 to a few times 
10 5 for 1 to 10 MeV ions, independent of the density of 
the detector medium. 

When multiple collisions are taken into account, the 
size of the neutron tracker for efficient detection to ac- 
commodate n collisions can be estimated as the following. 



R(n) = £\ 



(46) 



with R(n) being the total range and Sj the neutron range 
prior to the (j + l)th collision. Following [16], 



(s )(cos9i) 



(47) 
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FIG. 6. Required spatial dynamic range £> as a function of ion 
energy. Here 51 = 5L s /y/2 with 5L S being longitudinal range 
straggling. 1 H is in the methane gas. 2 H is in the deuterium 
gas. The results are not sensitive to pressure. 



with (cosGi) = 2/ (3 A). Therefore 



(R(n)) = (s ) ]T 

3=0 



3.4 



(48) 



Drop the vector notation and use S(n) for the aver- 
age detector size for n collisions, from the magnitude of 
(sn) = An, we found 



S(n) < S(oo) 



An 



1 - --' 

1 3.4 



(49) 



For 1 H, S(2) = 5A /3, S(3) = 19A /9. For 2 H, S(2) = 
4A /3, S(3) = 13A /9. For 4 He, S(2) = 7A /6, S(3) = 
43Ao/36. The deviation from the average detector size 
can be defined as 



ai(n) = (R 2 (n))-S 2 (n). 



(50) 



We are only interested in n < 3, which can be calculated 
readily. 



<4(1) = *o, 



a 2 (2) = A 2 3 



and 



4(3) 



An I - 



4 
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4 
9A 2 



9A 2 ' 



16 
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(51) 



(52) 



(53) 



IV. TRACK MEASUREMENT 



There are at least two types of ion trackers that are use- 
ful to real-time tracking of fast neutrons. One is based 
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on charge collection, as in various multi-wire proportional 
chambers (MWPCs) [17-19]. The other is based on pho- 
ton collection, as in optical imaging. The purposes are 
the same: a.) to pinpoint the positions of the neutron 
clastic collisions or locate the vertices A, B, and C in 
Fig. 1 and b.) to find the lengths and the directions of 
the recoil ion tracks, Qi and Q2- The number of charged 
pairs and photons produced during the stopping of the 
recoil ions determines the signals available. Both the 
charge and photon methods require single-particle sensi- 
tivity, that is, sensitivity to a single electron or a single 
photon. 

A. Charge- based ion tracking 

MWPCs deliver superb charge sensitivity and spatial 
resolution [17-19]. The latest high gain operations can 
achieve a single-stage charge multiplication by a factor 
of 10 3 to 10 6 , sufficient for single charge detection using 
low-noise preamplifiers. Fine spatial resolutions less than 
100 /jm are possible with micro-pattern readouts, such as 
micromegas and gas electron multipliers (GEMs) [19]. 

In the absence of a magnetic field, the spatial resolution 
in an MWPC (aa) is limited by the diffusion of primary 
electrons [20], 

/4e e £ 

ffd = fev (54) 

here Eb is the bias electric field strength in the drift re- 
gion, £ is the drift distance, e e is the so-called characteris- 
tic energy of drift electrons. Adding a magnetic field may 
reduce diffusion across the field lines as in time projection 
chambers [21, 22]. However, besides the overburden of 
bulky coils or heavy magnets, the benefits of a magnetic 
field are further compromised here since the uncertain- 
ties of ion tracking also depend on the position accuracy 
along the magnetic field, which is not improved by the 
magnetic field. In the thermal limits, e e = 3/2ksT e . In 
another limit, when e e — 3gdXeEb/4: ^> ksT e , one has 
®d = V.9dA£ with a coefficient gd on the order of one. 
Typical resolution ranges from less than 100 fim to sev- 
eral hundred /zm, depending on the drift distance £. The 
maximum drift distances for different 1 MeV recoil ions 
have been listed in Table. I. It is found that the drift 
distances are on the order of 10 cm under the standard 
temperature and a pressure at 10 bar. Charge diffusion 
will make it impractical to recover ion tracks for longer 
drift distances. 

The spatial dynamical range D is about 500 for a 10- 
cm 10-bar detector cell. A large cubic tracker with at 
least 20 cells each side would be required to reach a 25 of 
10 4 . Close to 100% intrinsic neutron tracking efficiency 
requires 20 3 <~ 10 4 cells. For fewer cells, the detection 
efficiency decreases and is proportional to the number 
of cells. A detection efficiency of 0.1% will only need 
tens of detector cells. For a detector with a few cells, the 
intrinsic efficiency is in the range of 10 -6 to 10 -4 . In this 



case, to pinpoint a neutron source requires an increase in 
detection time to compensate for the reduced intrinsic 
efficiency. 

TABLE I. Estimates of the maximum drift distances before 
the initial charge cloud becomes too diffusive. No magnetic 
field is assumed. The number of the elemental charge after 
drift is 10 in a square pixel with a size l p . A is the recoil ion 
range at an energy of 1 MeV. 5L S is the range straggling. W 
is the energy required to produce each charge pair. All gases 
are assumed to be at 10 bar STP. 

Recoil ion a 1 H 4 He 

X (mm) 2.36 5.88 2.84 

SL S (mm) 0.10 0.24 0.10 

V b (V) b 1000 1000 1000 

W (eV) c 29.1 36.4 46.0 

l p (fim) 200 200 200 

£ (cm) 10 9.1 8.1 

a X H in the CH 4 gas. 
b V b = E b Z. 

c Data from [20]. W's for 2 H and X H are assumed to be the same. 

An MWPC normally only measures the projection of 
a 3D ion track onto a 2D plane. To obtain the third 
coordinate along the drift direction, a time projection 
technique may be used. Time projection measurement 
based on the light emitted from the ion tracks will be 
discussed in Sec. IV C. 



B. Light-field imaging 

Ion tracking based on light emitted along the ion tracks 
has yet to be demonstrated [2] . Besides electron- ion pair 
creation along its track, a recoil ion can also induce light 
emissions by exciting electrons from ground states. De- 
excitations emit light. Another type of light emission is 
due to energetic electrons produced by the ion. Head- 
on collisions between the ion and electrons produce the 
most energetic electrons called <5-rays by J. J. Thomson. 
Lower energy electrons are ejected when the ion and elec- 
trons interact farther away. Except for <5-ray-induced 
light emissions, which can be mm or farther away from 
an ion track, most light emissions come from lower en- 
ergy electrons, which closely follow the ion track and can 
be used for ion tracking. Optical methods remove track 
blur due to charge diffusion. 

The problems, sharing the same flavors as in charge- 
based methods, are to determine the positions and 
lengths of the light-emitting columns due to ion stopping. 
Although the fine spatial structures of the light columns 
smaller than ion range straggling is not essential, the lim- 
ited number of available photons is a challenge. In-situ 
photon multiplication is not considered. We also neglect 
light scattering and absorption before reaching the de- 
tectors. 

Light-field imaging (LFI) extracts scene information 
from different properties of a photon, its direction, wave- 
length, arrival time and polarization [23, 24]. LFI can 
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minimize the number of photons required to localize a 
photon-starved scene. As few as two crossing views are 
sufficient to locate a point source [25-27]. Any addi- 
tional detail from the scene, such as the source structure, 
requires more photons [28]. Photon energies or wave- 
lengths dictate detector material use and construction for 
high detection efficiency. Photon arrivals at the detectors 
are randomly delayed from the times of excitation. The 
random photon emissions follow the Poisson distribution 
with time constants equal to the inverse transition prob- 
abilities, or the half lives of the excited states. 

We may estimate ideally, how many photons are 
needed to determine an ion track, which is approximated 
by a line segment, shown as the shaded region in Fig. 7. 
The question is to locate the two end-points. Equiva- 
lently, one may also locate an end point and measure the 
track length and the direction, (xq, yo, zq, sq, to, 1). Here 
the first three parameters are the coordinates of a point 
in the ion track, (so, to, 1) is the direction of the ion 
track. Each photon, represented by a ray as (u\, v\, w\, 
at, ti, 1), has to have an anchor point (xt, y\, z\) in the 
ion track. Each photon ray therefore supplies four inde- 
pendent equations that describe the fact that the anchor 
point has to be on the photon ray and in the ion track si- 
multaneously, Eqs. (55) and (56), in which the unknowns 
are (x , y Q , z , s , t ) and (xt, yi, Zt). 



(u1,v1,w1,s1,t1,1) 




(x2, y2, z2) 

(x1, y1, z1 



(U2,v2,w2,s2,t2,1) 

FIG. 7. An ion track is approximated by a line segment 
(shaded). Every photon emitted, represented by a ray, has 
to anchor in the ion track. In general, the rays are not copla- 
nar. 



and 



Xt ~ Xq 

so 



X\ — U\ 
Si 



yi - y o 

to 



yi -vi 
u 



= zi 



zo, 



Z\ - Wt- 



(55) 



(56) 



Another four photon rays, or for a total of five photon 
rays, would have a sufficient number of equations to solve 
for a total of five unknown anchor points and the un- 
knowns for the ion track. For the track length, one needs 
to trace an additional photon. 

Although less than ten photons are sufficient for ion 
tracking, it will take more photons to locate the end- 
points of the ion track and measure the ion track length 



more accurately. The least number of photons needed is 
determined by Poisson statistics. For a signal-to-noise of 
five, the minimum number of photons needed per resolv- 
able segment of the ion track is 25. If the whole ion track 
is considered, the total number of photons is to multiply 
25 by the number of resolvable segments, which is the 
ratio of the ion track length to straggling length, or 25 
From Fig. 5. Therefore the total number of photons is at 
least 625 for MeV recoil ion tracks. 

Photons emitted along an ion track come from 
two sources of electrons. Multiply scattered electrons 
(MSEs) [29-32] and 6-rays. MSEs have lower energies 
than (5-rays and form the primary track [33, 34]. The 
relative numbers of MSEs to 6 rays is proportional to the 
ratio of multiple scattering cross section to two-body col- 
lision cross section. The ratio may be characterized by a 
A parameter [35, 36], 

A S £. (57) 
Here, the so-called shielding length (Xd) is given by 



A; 



noZe 2 



(58) 



with Io being the ionization potential of the stopping 
medium, no the atomic density of the medium, Z the 
ionic charge and eo the vacuum permittivity. The dis- 
tance of the closest encounter b is given by 



b 



Ze 2 



AireoE, 



max 

<5 



(59) 



^max 



The maximum energy of a S-iay is given by Eg 
4(m e /M p )(Ei/Ai), which is the largest for protons with 
Ai = 1. At the present, E l < 10 MeV. Therefore, 
E™ ax < 22 keV. Both ranges of MSEs and 5 rays may be 
calculated from the continuous slowing-down approxima- 
tion (CSDA). The lateral distances of S ray photons from 
the ion trajectory is larger than MSE photons. In the hy- 
drogen gas, the CSDA range for a 22 keV electron is 4.54 
xlO -4 g/cm 2 [37], or about 5 mm at 10 bar pressure and 
room temperature. 

The primary track has a sharp radius, which is de- 
termined by MSEs with an average energy up to 1 to 
4 keV [33]. The primary column radius, the same as 
the average MSE range, may be approximated by R e = 
9.93 x 10 -6 i? e [g-cm _2 keV _1 ], independent of stopping 
medium [38]. 

The number of photons per unit length along an ion 
track is given by 



dNx 
dL 



= -(l-F) 



1 dE f((Jiv e 

W dL (<TiVe) 



(60) 



at a wavelength A. W is again the average energy for 
ion-pair production, which is assumed to be the same 
for both MSEs and 5-rays. (<Jiv e ) / ((JiV e ) is the ratio of 
the averaged excitation rate to an atomic level / to the 
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ionization rate. / < 1 is the branching ratio for A. F 
is the fraction of energy that goes to the <5-ray produc- 
tion and does not contribute to the light production in 
the primary track. F was found to depend on Ei/Ai 
only [33] . After an ion energy falls below a certain thresh- 
old, (5-rays can no longer be produced, and F = 0. From 
Ef" = A(m e /M p )(E i /A l ) < 1 keV, ones finds E t = 
0.46 MeV for protons, 0.92 MeV for deuterons and 1.8 
MeV for 4 He's. These numbers are small compared with 
the values in Fig. 3 of [33] . If one uses the threshold of 4 
kcV instead, then E { = 1.8 MeV for protons, 3.6 MeV for 
deuterons and 7.2 MeV for 4 He's. Here we use F = 0.25 
for simplicity, which may underestimate the photon pro- 
duction by up to 25%. Furthermore, Eq. (60) does not 
include photons from ion-pair recombination. 

Based on the references [39-41], the branching ratios 
are calculated for several wavelengths of 1 H, 2 H and 4 Hc, 
Table. II. The spectral properties of 2 H and 1 H are 
treated the same and not repeated. The electron en- 
ergy distribution in the primary track is assumed to be a 
gaussian with an adjustable temperature parameter. The 
resulting ratios (<Jiv e ) / (<TiV e ) are shown to be insensitive 
to the electron temperature, Fig. 8. 



TABLE II. Calculated branching ratios (/) for several wave- 
lengths of the 1 H and 4 He gas. The largest numbers of pho- 
tons are expected at these wavelengths. The results for the 
2 H gas are expected to be the same as those of X H and not 
listed. 



Wavelength (nm) 


Transition 




;io 8 s- 1 




H I 121.6 (L a ) 


2p ^£3/2,1/2 


Is 


S'l/2 


4.70 


1.0 


H I 656.3 (H a ) 


3p ^3/2,1/2 


-> 2s 


S'l/2 


0.441 


0.44 


H I 102.6 {Lp) 


3p £3/2,1/2 


Is 


Sl/2 


0.558 


0.56 


He I 58.4 


2pls 1 P 1 - 


► Is 2 


^0 


18.0 


1.0 


He I 501.6 


3pls 1 P 1 -> 


2sls 


"So 


0.134 


0.023 


He I 53.7 


3pls 1 P 1 - 


► Is 2 


'So 


5.66 


0.98 



a Transition probability. 

b Collisional depopulation rate is neglected. 

Combining the results from Table. II and Fig. 8, we 
obtain the expected photon densities along the primary 
tracks, Fig. 9 for X H ( 2 H) and Fig. 10 for 4 He. Spectral 
lines H I 121.6 nm, H I 656.3 nm and He I 58.4 nm have 
sufficient photons per straggling length for imaging. He I 
501.6 nm is not a good candidate because of the paucity 
of photons. Using lense-based method to measure pho- 
ton directions (see below) , we will also run into problems 
with He I 58.4 nm, since the shortest transmission wave- 
length of any known material is about 105 nm using a 
LiF window or lense. It is worthwhile to look into lense- 
less methods [42] to detect directions of VUV photons 
like He I 58.4 nm, but we do not elaborate here. 

A simple 47r detector without any optics can capture 
all the light but is not useful because it can not resolve 
the directions of photons. On the other hand, a detector, 
in conjunction with small pin-hole collimators, can meet 
the angular resolution requirements but has too small 



1.2 



I 1 '° 
to 

N 

c 0.8 
o 

■| 0.6 

ro 

o 

i« 0.4 



0.2 
0.0 



H I 121.6 nm 



He I 58.4 nm 



H lfiS B3nm, HI mpfinm 

He I 501 6 nm, Ha I 53 7 nm 



10 20 30 40 50 

Electron temperature (x13.6 eV) 

FIG. 8. Calculated ratio (criv e ) / {<7j« e ) for a gaussian energy 
distribution of electrons in the primary track. The results are 
insensitive to the electron temperature. We use 0.91 for H I 
121.6 nm, 0.24 for He I 58.4 nm, 0.15 for H I 656.3 nm and H 
I 102.6 nm, 0.058 for He I 501.6 nm and He I 53.7 nm. 




1 .0 

Ion energy (MeV) 



10.0 



FIG. 9. Calculated photon counts per lateral straggling length 
in the primary track for 1 ¥L and 2 H. 1 ¥L is in pure hydrogen 
gas. 2 H is in deuterated methane. The results are not sensi- 
tive to the mass density (pressure). 



a photon collection efficiency to be useful. A possible 
method is to use a lenslet array that covers ~ 4-7T solid 
angle so that most photons are captured. For each lenslet 
in the array, there is a two-dimensional detector array 
at the focal plane, Fig. 11. Each lenslet in the array 
sorts the incoming photons into different detector pixels 
according to photon directions. The pixel size of the 
detector array (l p ) matches the diffraction limited spatial 
resolution, l p = a r f p , with f p being the focal length and 
the angular resolution a r = 1.22X/D [43, 44]. Here D is 
the diameter of the lenslet, as in Fig. 11. 

The spatial resolution at the scene is given by 

61 = a r O + D = 1.22^— + D, (61) 
with O being the distance between the ion track and a 



10 



10" 




He I 501.6 nm 



1.0 

Ion energy (MeV) 



10.0 



FIG. 10. Calculated photon counts per lateral straggling 
length in the primary track for 4 He. 4 He is in pure helium 
gas. The results are not sensitive to the mass density (pres- 
sure). Only the He I 58.4 nm has sufficient photons above the 
threshold of 10 counts per straggling length. The visible line 
He I 501.6 nm has too few number of photons. 



C. Time measurement 

Recording time information is useful and in some cases 
necessary in fast neutron tracking. As in Sec. II, the num- 
ber of required recoil collisions reduces from three to two 
if neutron TOF is available. Timing of neutron collisions 
can be used in charge-based ion tracking to obtain the 
ion drift time and the third coordinate along the drift 
electric field. In LFI, photon TOF from an ion track 
to the detectors is redundant because of Eqs. (55) and 
(56). We shall only discuss timing measurements based 
on light emissions. Other methods, such as sensing neu- 
tron magnetic momentum, may also be possible but are 
not discussed. 

For a photon emission (A) with a transition probability 
Ai\, the number of atoms in the corresponding excited 
state (Ni) as a function of time is given by 



dNjjt) 
dt 



A' 



(1-F) 



1 dE (aiv e 
W dt (<JiV e 



(64) 




Summing over different wavelengths is due to the fact 
that several de-excitation transitions may be possible. 
The number of photons with the wavelength A as a func- 
tion of time N\(i) satisfies 



dN x {t) 
dt 



AaN^t). 



(65) 



Here t = is when a neutron collides with the ion, 
which coincides with the start of ion motion. Combining 
Eqs. (64) and (65), 



FIG. 11. An approach to photon-starved LFI is to use a lenslet 
array that surrounds the scene, so that close to 100% photons 
are detected. Each lenslet sorts photons into different pixels 
of a detector according to their incoming directions. The 
directional uncertainty is diffraction limited. The possibility 
of super-resolution is not considered. 



dNijt) 
dt 



ldNx(t) {1 _ F) J_dE(^!A. (66) 



/ dt 



W dt (<7iV e 



Eq. (66) is consistent with Eq. (60) by integration over 
time from the beginning of the ion motion until ion stop- 
ping or the end of de-excitation time. Solving Eq. (64) 
fior JVi(t), 



lenslet. The best resolution of 5l min — 2D is therefore 
reached when 



-V,(/)= / dse^-Hl-F)±-fp^, (67) 

W ds (a l v e ) 



D 1 = 1.22AO. 



(62) here we define 7j = J2\Ai\- Solving Eq. (66) for N\(t), 



For Sl min = 5L s /\/2, 



D = 



SL S 
2V2' 



(63) 



For 5L S — 200 /zm, D =71 fim. The distance O is in- 
versely proportional to the wavelength. At H I 121.6 nm, 
O = 34 mm. At H I 656.3 nm, O = 6.2 mm. At He I 58.4 
nm, O = 70 mm. Correspondingly, the spatial dynamic 
range of a detector cell, D = 20/5L s , is 337, 62 and 702. 
On average, no more than one photon will reach each 
lenslet. 



Nm = m-F) m - m M -m(t). (68) 



w 



{cTiVe} 



The slowing-down times of recoil ions are comparable or 
longer than de-excitation times, which are around 1 ns, 
see Fig. 12. Here the slowing-down times (tsl) are given 
by 



TSL 



E(0) 



dE f dE 



(3c \dL 



(69) 



where dE/dL is the energy loss of the recoil ions per unit 
length, j3 is the ratio of an ion speed to the speed of light. 
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FIG. 12. Average ion slowing-down time in different gases 
as a function of the initial ion energy. From the top to the 
bottom are for electrons in 4 He, 1 H in 1 H gas, 2 H in 2 H gas, 
1 H in methane gas, and 4 He in 4 He gas. All gas pressures are 
at 10 bar STP. 



FIG. 13. Ion energy loss rate in different gases as a function 
of the initial ion energy. From the top to the bottom are for 
4 He in 4 He gas, 1 H in methane gas, 2 H in 2 H gas, and 1 H in 
1 H gas. All gas pressures are at 10 bar STP. 



We used SRIM to calculate dE/dL for different ions [14] 
and ESTAR database maintained by NIST to calculate 
dE/dL for fast electrons [37]. 

The time of neutron collision (t=0) can not be mea- 
sured since there is no signal available. The emission of a 
photon obeys Poisson statistics, so there is a theoretical 
possibility that at t — a photon with the wavelength 
A is emitted and is detected after a time delay which is 
the sum of photon TOF and detector response time. The 
average time that characterizes the Poisson statistics of 
photon emission to satisfies 



TABLE III. Calculated to for various scenarios. The ions are 
at 2 MeV. 



Wavelength (nm) 


7i (ns x ) 


T ex (ps) t (ps) 


H I 121.6 {L a ) 


0.47 


0.03 


11 


H I 656.3 {H a ) 


0.10 


0.18 


92 


H I 102.6 {Lp) 


0.10 


0.18 


81 


D I 121.6 


0.47 


0.09 


20 


D I 656.3 


0.10 


0.54 


157 


D I 102.6 


0.10 


0.54 


139 


He I 58.4 


1.8 


0.79 


9.4 


He I 501.6 


0.58 


0.33 


222 


He I 53.7 


0.58 


0.33 


34 



N x (i ) = 1 



(70) 



with N\(t) given by Eq. (68). Using the approximation 
dE/dt =const, which is justified a posteriori , 



to = 



with r ex being defined as 



2r 

^ ' ex 
fit 



(71) 



W (a lVe ) (dEX 



1 - F (atv e ) \dt ) 



(72) 



Using data from Fig. 13 and Table. II, we may cal- 
culate to for different scenarios. The results are given 
in Table. III. To take advantage of the ps time re- 
sponse of the photons, efficient, single-photon sensitive 
and ps-resolution detector would be ideal. These de- 
tectors do exist and their performances are suitable for 
LFI [45, 46]. If 100 ps time resolution is indeed obtained, 
from Eq. (44), the neutron TOF for 8% resolution will be 
1.25 ns. At 10 MeV, the flight length required is about 5 



V. DISCUSSION 

Our analysis has mostly focused on elastic collisions 
of MeV neutrons with 1 H, 2 H and 4 He, all of which are 
abundant and can exist in the gas phase. For *H and 
2 H, abundant organic compounds also exist. There are 
a few low-Z elements that have relatively large absorp- 
tion cross sections for MeV neutrons. In Fig. 14, the 
ratios of the fast neutron absorption cross sections to the 
elastic cross section of 1 H(n,n) 1 H are shown for 3 He, 6 Li 
and 10 B. 3 He is the most attractive in that it has the 
largest absorption cross section and the smallest Q-value 
at 0.76 MeV. However, due to the short supply of 3 Hc, 
a large volume of 3 He may not be available. In a small 
volume, with the detector areal density (nL) satisfying 
anL <C 1 and a being either the neutron absorption or 
elastic collision cross section, the efficiency of an absorp- 
tive neutron tracker can be higher than that of an elastic 
tracker. Since an absorptive tracker only requires one 
collision, the intrinsic efficiency is given by t\ ~ o~ a b s nL. 
For an elastic tracker that requires three collisions, the 
intrinsic efficiency €3 ~ (a e inL) 3 . For an elastic tracker 
that is based on two collisions, the intrinsic efficiency 
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£2 ~ (o'einL) 2 . MWPCs and LFI can be used as absorp- 
tive trackers with 3 He, 6 Li and 10 B. Some other atoms 
have to added to a gas mixture so that photons with ap- 
propriate wavelengths and sufficient intensities are avail- 
able for LFI and time measurements. Atoms or com- 
pounds of 3 He, 6 Li, 10 B, together with a sources such as 
241 Am, will be useful to developing LFI technique. 

0.4 P 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ^1 



0.3 




o.o L_ , , i , , , i , , , i , , , i , , 

2 4 6 8 10 

Neutron energy (MeV) 



FIG. 14. Ratios of the neutron absorption cross sections of 
3 He, 6 Li and 10 B to the elastic cross section of 1 H(n,n) 1 H for 
neutrons in the MeV energy range. The data are from [13]. 

Detector pressures here are in the range of 1 to 10 
bar STP. Although higher pressures or even solid density 
trackers may be desirable for higher intrinsic efficiency, 
higher densities may not be compatible with the MWPC 
operation. In LFI, solid density would require optics 
with sub-/zm spatial resolution, which may compound 
the technology development in the beginning. Another 
difficulty is the apparent lack of lense material to focus 
VUV photons such as He I 58.4 nm. Initial experiments 
with 1 H or 2 H is easier, when both the visible line H I 
656.3 nm and VUV line H 1 121.6 nm are available. Other 
wavelengths may also be possible by using different gas 
mixtures, although we have only discussed H2, CH4 and 
their deuterated versions here. 

Spatial resolution is limited by ion range straggling. 
For individual ions, however, the ultimate spatial resolu- 
tion of ion tracking is limited to the primary ion track 
width, or the ranges of multiply scattered electrons as 
given by R e = 9.93 x 10~ 6 E e [g-cm _2 keV _1 ], indepen- 
dent of stopping medium. For electron energies E e = 1 to 
4 keV, R e = 1 to 4 mm in 1 bar STP hydrogen gas (mass 



density 8.93 x 10~ 5 g/cm -3 ). At 10 bar pressure, the 
range shrinks by a factor of 10. In other gas mixtures, 
such as CH4, the primary ion track width will be smaller 
by a factor of 7 from that of hydrogen. In all cases, the 
ion track widths are comparable or less than the range 
straggling of ions with energy above 1 MeV. 

One possible application of the fast neutron trackers is 
to pin-point the locations of unknown neutron sources. 
Compared with a recoil proton telescope [10, 47], the 
proposed tracks have better angular resolution. Further- 
more, when the TOF of neutrons can be measured and 
therefore two recoil events are sufficient for each fast neu- 
tron, the intrinsic efficiencies of the proposed trackers will 
be comparable to a recoil proton telescope with a similar 
areal density. A significant reduction of neutron events 
may be possible in order to pin-point an unknown neu- 
tron source using the proposed trackers. 



VI. CONCLUSIONS 

We have discussed the principle of fast neutron track- 
ing based on elastic collisions. The linear momentum 
of a fast neutron can be measured from as few as two 
and a 'half consecutive recoil ion tracks. If the time- 
delay between two consecutive recoil ion tracks is also 
measured, the number of ion track measurements can be 
reduced to one and a half. The magnitude and angular 
resolution of neutron momentum are limited by ion range 
straggling and ultimately by ion track width in the stop- 
ping medium. The angular resolution and magnitude of 
neutron moment measurement are about ten percent due 
to ion range straggling. Multi-wire proportional cham- 
bers and light-field imaging are proposed for fast neu- 
tron tracking. Single-charge or single-photon detection 
sensitivity would be required. Light-field imaging is free 
of charge-diffusion-induced image blur, but the limited 
number of photons available can be a challenge. 1 H, 2 H 
and 3 He would be used for the initial demonstration of 
fast neutron trackers based on light-field imaging. 
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